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THE DEVELOPMENT OP THE POINCARE-SIMILAR ELEfiENTS WITH 
THE TRUE ANOMALY AS THE INDEPENDENT VARIABLE 
By Alan C. Mueller 

Analytical and Computational Mathematics, Inc. 


1.0 INTRODUCTION 


In reference 1, Scheifele established the Hamiltonian of the unperturbed 
two-body problem in extended phase spacer Dependinj? on the type of time 
transformation, eight canonical elements were develoned with the true 
anomaly or the eccentric anomaly as the independent variable. These two 
new sets, DS(*) and DS(u), however, contain singularities for small eccen- 
tricities and inclinations. In reference 2, these singularities are removed 
by a transformation from DS(u) to eight canonical PS(u) elements. As 
in the DS(u) variables, the PS(u) variables have the eccentric anomaly 
as the independent variable . 

In reference 3» tbe DS(*) variables are transformed to the PS(*) 
elements to remove the singularities. However, no direct relation was 
established between the eight canonical PS(*) elements and the Cartesian 
coordinates. It is the purpose of this report to establish those relations 
and to develop the perturbed equations of motion in the PSC*^) space. 

As will be seen, the relationships are not trivial to establish; however, 
once obtained, they are found bo be rather simple expressions.' 

Lastly, this report will demonstrate the accuracy of this new set 
when it is applied to numerical oi>bit prediction problems. 
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2.0 TRANSFORMATION FROM DS(* ) TO PS(* ) VARIABLES 


In reference 1, Soheifele presents the Delavnay-simllar (DS) elements 
with the true anomaly as the independent variable - DS(*). The angular 
variables are as follows. 

^ s the true anomaly. 

g = the argument of perifocus. 

h = the argument cf >-he ascending node, 

i = the time element. 

The action variables are as follows. 

♦ = conjugate to ♦ , related to two~body energy. 

G = the {iotal angular momentum. 

H s the Z-component of the angular momentum. 

L = the total energy (two-body energy plus perturbing potential). 

The corresponding Hamiltonian with a perturbing potential V is given 
by 

( 1 ) 

( 2 ) 


F s ❖ - 




r2 

— V 


where 


° " 2 * 


The equations of motion can be written as 


'' d4 a F 

— = f = *, g, h, 1 

dt aL' 


dLt 3F 

d— = -nT H, L 


Note that in unperturbed motion (V = 0) 

l|> = T + * _ 


(H) 


i = 


(2L„) 3/2 


X + I, 


where t is the independent variable. All other variables remain constant. 

The DS( «)-elements , however, become ill defined for vanishing eccentric 
ities and inclinations. To remove these singularities, one should transform 


2 


I 'i 


to a type of element aJmilar to the Poinoane variables. These ?S{^) elements 
will bo named o (coordinates) and js (momenta). As in the DS(«) 
elements, the P'S(V) elements have the~true anomaly a:^ the dependent variable. 
Those elements should not be confused with the PS(u) 'variables oresented 
in reference 2 ; howovor, In many aspects they are similar^ 


The PS(^ ) elements can be defined by a canonical tbansformation from 
the DS(* ) elements by the generating function S. 


1 O 

t p j| + (♦ + g + h)p 1 - 5 P2_Jfcan(g + h)] 


- ^ p 3 ^ tan (h) 


where 


as , . as 

OB — , L' = — - 

“ , ' a I ' 

Dp 


and 


t? _ 

,T . 


0 * = (oi, 0 2 , 03 , Olj) 


_P ^ = (p 1 , 0 2 • ** 3 * ^ ^ 

The transformation yields the coordinates, 
0|-^ -fK + h, 

P 2 = - \j2{» ~ G) sin(g + h) 

03 = - \j2{G - H) sin(h) 

0^ = t 

and the momentas, 

_ . if 

P -J - ^ 


P 2 =\(2(« - G) cos(g + h) 


P 3 =\/2(G - H) oos(h) 

P i| = L 

The Hamiltonian now written in PS(<|>) elements reads 

,2 


F 


Pi - ^ + 
>/2o ii 


(5} 


( 6 ) 


(7) 


( 8 ) 






3 








3.0 EQUATIONS OF MOTION 


In roference , a method la presented in which the equations of motion 
in extended phase space are extended to include forces that cannot be 
defined as derivatives of the potential. The equations. of motions in 
Cartesian space may bo written as 




£ 


n (12) 


where F is the perturbing force, V is the perturbing potential, and 
£ is the force not expressed as a potential. If the Unperturbed Hamiltonian 
la expressed as 


^0 ■ 1 




(13) 


then the equabiona of motion in PSI^) space now become 


do ^ 
dr ip 


1= r + V ^ ^ ^ (il\ 


(HI) 


^ = V + 1- * s 

Ot 3^ q afl “ 


r^ av 


where S’ aad £ are the canonical forces and are defined as follows. 

B = d at p A , u 
q 3p 


5 = ii: pT p , B = 
q - ’ 3 0 


(15) 


aV 


aV 


a p ax 


= bT 


-■? 3 ff 


ax 
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P‘ 


P4 s - + X2*P2 + X3*P3) 

Xg, X3 = velocity vector 
Equation can now be written aa 



The roatrioea A and B can be found from the equations relating the 
Cartesian coordinates and time in terms of 0 and 0, Here, also, it 
is worthy to note that i “ ~ ' 



is the time transformation equation. 
‘:o Integrate to find the time. 


However, this equation is not ua. 



il.o CARTESI/jircOOnDINATBS (xi, Xg, X3, t, x,, Xg, X3) 
IN TERMS OF P5(*) ELEMENTS (e, 0) 


Although the PS(*) elements are uniquely defined by equations (6) 
and (7), their relations with the Cartesian state vector must be given 
directiy, without making use of the DS(^) eloinents, since ^ and g are 
numeric^aily ill defined for near-oircular orbits and h is ill defined 
for small inclination orbits. However, the DS(^) elements G and H 
may be used occasionally as abbreviations since they are always well defined. 


Prom reference 1,, the coordinates and the time can be obtained from 
the DS{^) elements as 


where 


x^ = r Ccos(* + g)co3 h - sin($ + g)sin h cos I] 
Xg,,= c [oos{* + g)sin h + sin(* + g)cos h cos I] 
X 3 = r sin(i + g)sln I 

u 


Xj| = t = j, + 


(E - * - ” yr~-T^ e sin ♦) 

p V 


07) 


r = 


1 + e cos 


= :(“-•* ar) 


7T 


2L 


08) 

(19) 

( 20 ) 


C03 I = g 


E = 2 arc tan 


(fi 


- a tan t 
e 2i 


( 21 ) 

( 22 ) 


'Note here that these notations are used as abbreviations in the DS(>) 
theory to transform from the coordinates to DS(,t) elements. Although 
they are similar in form to Keplerian elements, they are the same only 
when the perturbing potential V = 0. 

By trigonometric Identities, equations 07) may be rewritten as 
x^ = r sin^ g cos [h - (* + g)] + r cos^ - cos(^ + g + h) 
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T X 

xg = r sin^ i sin [h - (* + g)] + n oos^ | sln( * + g + h) (23) 
I I 

Xg s 2 r cos 5 sin = sin (* + «) 



From equations (5) and ( 6 )» it is oasiiv shown that 
* - Q B - (02^ 

and 

G - H s I ^^3^ ■*■ *^3^^ 

Then » using equations (21) and (2^1), we find 
\/ 2 (G - H) s 2 ^/^ sin I = + 03 ^ 


Equation 19 may be wnittr^-' >s 


P vi f'" 1 ^•’2“' + “2^^ tI; 


2 


From the identity 

sin o| sin (g + h) + cos cos(g + h) = 00s i 
We find by mulbiolication with ^z(i ^ G) 

- og sin O'! + eg I cos 0^ = \jz(^ - G) cos ^ 

Similarly, one can derive 

eg sin 0^ + og cos 0^ = \J^V “ G) sin 
From the identity 

sin oi sin h + cos 0^ cos h = cos (^ + g) 
one can derive 

-03 sin 01 + P3 cos 0^ = 2 v/G* sin | cos(* + g) 
and, similarly, 

03 sinvo^ + *>3 Gos oi = 2 v/G sin g sin(* + g) 

8 



(2H) 


( 25 ) 


(26) 


(27) 

(28) 



(30) 



Although ^ is ill defined for small eccentricities, the right-hanri sides 
of equations (27) and (28) remain well defined because multiplication 

factor* v/2(« - g) becomeO quite small. The same argument applies’’ for 
equations (28) and 29) if the inclination is small where |(> + g)) is 
ill defined. 


By using equations (20), (27), and (28), one can derive expressions 
e cos = 7.iQ '(31) 

(32) 


where 


and 


e sin = ZgO 

Zi = P2 cos 01 - og sih 01 
Zg = 02 ®2 


1/2 


(33) 

(3*0 

(35) 


There is no numerical problem in computing 0 except in the case where 


e 1 or 


2 u ..1 / 2 ^. 2 - 




va 5 (o2 + P2 > 


• ^ We can now make use of equations ( 6 ), (7)t (29)* and (30) to substitute 
in equations ( 23 )* - 


x*| ^ 1 

^2 ” ^3 *** ^ sin O') 

X3 = R* /ITgTIo 


(36) 


where 


B* = g 

2G 


and 


R s" p 2 s in 0 1 +03 cos 0 1 


( 37 ) 


(? 8 ) 


Note that the second terms in Xi and Xg are dominant for small eccentri- 
cities and inclinations. 


The redius, 


(39) 


P 

r 5 — 

-1 + e cos ^ 

-is evaluated by using equations (26) and (3D* 
The time equation, 


Xh = t = oil + '^"TTo (E - ^ ^ )/l - e sin »> (1<0) 

(2p4)3'2 D, 

is evaluated by using equations (20) and (32) and by - 

E - <> = -2 arc ban / ^ A (*» 1 ) 

y 1 + \/1 - + e cos ^ / 

We may now proceed to the comnutation of the velocity components. 

In the case of the DS(<^ ) elements, the velocity components were preferably 
computed by forming the time derivatives of the transformation equations 
for the coordinates. The reason was bo avoid the numerical difficulties 
that are encountered at some pfiinbs of the orbit. These numerical dif- 
ficulties were caused by the irtbermediate step of computing the spherical 
coordinates. 

The time derivatives of the 08(4 ) elements are different from zero 
and depend on the current value of the perturbing potential, Vg. We 
denote the "unperturbed orbit" that is obtained from the Hamiltonian, 


’^DS 



(« 2 ) 


where Vq is a constant, the "DS orbit." 

Because g, h, t, and H do not appear in the Hamiltonian, the 
DS('t') orbit has the following properties. 


G = const., H = const., L = const., h = const. ( 43 ) 

ThO means that the orbital plane and the total energy remain unchanged 
or I'iiat G, H, L, and h are osculating elements. By examination 
of equation (17), one must conclude that, in fact, the quantity 4 + g 
is osculating. Thus, from the angular momentum law, * : 


db " dt " p2 

If these facts are taken into account in the PS(4 ) element system, we 
can conclude from ecuations (6) and ( 7 ) that 



(•lH) 


( 45 ) 
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where (’)"'= ^ 


The time derivablvee of bho roinaiiiing HS(<>) olemenbs oan be oompubed 
by esbabliahing bhe canoniohl equablona bo bho Hamilbontan . 


U 

FpS = Pi “ “= + — V, 




.2 ^ 


s/2 Plj 


) 


I0„ 

02- + 02- - 01 - 
Firab , bhe equabions for 02 and 02 oababliahod, 


(>IG) 

(•17) 


l!B . !!p« -V 


d T Do 


0 3 0^ 


(48) 




dt ao/ 


= V 


(fj 


0 ao^ 


(49) 


Noting equation (lua), we find 

- V * 2 - 

- - - 0 ^,2 ao2 


^2 - i- 


db 


(50) 


g - ^°2,dt Q \*^/ 


2 - d;J db ~ Fo~2 


(51) 


However, the time derivabivea of 03 and Og are nob needed direobly 
bub aa oombiinatioha - as are 


Og 02+ Pg Pg 


(52) 


Before we proceed bo the computation of expression (52), we will observe 

(53) 


o p 

that if C is a function of bhe expression og^ + 03 “ 


b = b(og2 + og2) 


) 


and does nob denend on og or 03 in any other fashiop;t then 

= 0 -"V' 

for symmebrioal reasons «. 


3 b a c 

“2 17^ - “2 a7g 


(54) 


11 




i 


Observe that b. a, an<l Q are funetiona of the tyre found In enuatiun 
(53). 


Thus, 


. u ° 

0202 P2<’2 ''c ^2 



■(s' 


2 




3Q 2 

o o “ P 

^ 3P2 


Si)l 




Noting equations (32) and (39), we find 


2 r / 8P 9P_\ „ £I j, 

02°2 t «2“2 = '^c p B ( iTi ~ 302/ P ®’’2 


2 r a(e 003 *) 


- P 2 


3 (e cos iji ) 


= -V. r Z 


0 n ri “2 


+ Q fog 0®® 01 + P2 01 



s _|T Vq Q fog ®os 01 + pg sin 0 


Next, let us compute oi, where we note that oi occurs in r only y 
e cos * in equation (31). 

n sFdo fl 3 / r^\ “2 ,, 3r 


q V ~ ( — )= — Vq — 

M = ::? nr ■ .2 ® aoi Vo / r « ao 


' pdi ao 1 P'i «« 1 \'' / ' - ' 

2r a (e o®® » ) _ i£T Vq Q (pg sin o i + og cos o i j (56) 

p /'C 3o^ ^ \ > 

Finally, »e nrooeed to the oonnutation of the time derivative of the distance, 
f, on a "PS orbit.” For this hurpose. It la oonvenlent to Five the derivative 
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= — 0(^2 sin 0 ^ + ° 2 ° 1 ) 


To obtain q, we use eouations (55) and (56). 


/aq . 3q 
Q = 7 L — P 1, + r— 


a- 


3 O 1 


1 • . • * 

-p - 2 ° 2 ^ 2 ^^ 2 ) 


k =1 


- VpO (p2 sin a.] +03 ®°® ° 1) 


Finally, equation (58) is solved for r. 

2 ' 


. r . q r 
" 2q 2r I q 


and exoressions (59) and (60) are inserted, with the result that 




a O GOS o 1 + p o sin a 





Noting equations (31) and (33)f we find 




( 61 ) 


In DS theory (and, accordingly, 
vanish on the DS (or PS) orbit, 
for the Hamiltonian to vanish, 



PS theory) , one requires that the Hamiltonian 
Noting equation (46), one can see that 



( 62 ) 


But from equation (47), 


2q 



2(q - p^) + + pg2 = (63) 

Inserting equat-ions (;62) aiid ( 63 ) into equation (61), we obtain the final 
fi3i''ii for r. 


r 


e sin ♦ 
P - 




(64) 


We now can obtain the velocities by differentiatinp equations ( 36 ). 


• * i • * >1 

*1 “ - R "3 + r cos ®i - r®^ sii,\ ° ^ 
i <2 = - + r sin °i + r®-| cos ^ 

X 3 = R*V2(G + H) 



(65) 


where 



• • 

rR + Rr 


R = “i(f3 cos 0 ^ 


V 3 sin 











Noting-equation (*»5) and grouping terms, we arrive at the final form of 
the velocities. 


*1 


• n , G 

-R Oj + r cos - - sin °-| 


i! 




S* » • j n 0 D 

-R '*3 + r sin °i + - cos 


^3 = 


* II 

n ** 


R* y/2(G + H) , 

* • 

rR + Rr • 0 / 


(6 

( 6 ' 

(61 









5.0 PS(«) ELEMENTS ( ?, 7 ) IN TERMS OF THE ' 

CARTESIAN COORDINATES (x,, Xg, X 3 , t, x^, X 2 , X 3 ) 


From reference 3, we obtain 
Gi s X2X3 - X3X2 

G2 = X 3 X 1 - X 1 X 3 

G3 = Xix| - XgXi 

G = + Gg^ + G3' 

H = Go 



L !i 


2 l^1 


+ + xo2)- V 


where 


r = + xg2 + X3‘ 


Since the value of the value of the homoReneous Hamiltonian (eq.( 1 )) must 
vanish, one finds that, in DS((>) theory, 

♦ = G - v/g^ - Rr^V + -2i- ■■ (77) 

From equation (7), '' 

= <», Pjj = L (78) 

InsertinR these values in equations (2) and (19), one arrives at 

q = G - 5 ♦ + (79) 

2 am 


0 = J / G - « + ( 

And noting equation (2i|), we find that equation (35) can be expressed 
a s 


0 = 1 Pii 

u ^ 





■■ From expressions (31) and (39)» one obtains 


Zi = 


(P/r - 1) 


-1 - - Q 

and if we note that 


( 82 ) 







G P'1 - ^ (og^ + 0 ^) 

then one obtains from expressions (32) and (64) 


where 


^2 = Q(2a - g) 



(83)' 


(84) 


(85) 


By the definition of the angular momentum, its components can be Written 
as 

Gi = G sin 1 sin h - - 

Gg = -G sin I cos h (36) 

Gg = G oos I r 

From expressions (6), (7)i (25), and (86), one can derive 


= -- 

^ " C2(G + H)]1^2 
- 2G2 

[2(G + H)]^''2 

By eouations (36) and (37), we obtain 
2G(x3/r) 

" C2(G 
r cos 
r sin 

and thus 

ffi = arc 


= x-j -f R 0^ 
= Xg R 03 



(87) 

( 88 ) 

(89) 

( 90 ) 

(91) 


(92) 


Solvlnjj equations ( 33 ) and ( 3 ^J) for 03 and 03 1 we find 
03 = Z3 cos oj “ Zj sin (/ 

(>3 ~ Z3 sin O'! + cos ©I 


( 93 ) 

(9^0 


NobinK equations (20), (31) j (32), (*10), and (*I1), one obtains the final 
element i 


Oi} = t “ 


(20ij)3^2 


where 


E - ^ s “2 arc tan 


« - p.zaQ e^' 


Zs‘0 

( 1 + y/i^T^ + Z 1 • Q ) J 


) 


( 95 ) 


I ) 


(567 


and 


\/l - e^- = J— P 


M 


( 97 ) 


A sequence of the computations that sive the transformations to and from 
the elements and coordinates will be definitely outlined in the appendix. 











6.0 MATRICES A, B, 


> AND 




13 I3n.thy-b„t w?fi"?e‘5:u^?e^3““\SrJ:jSiJrare 


a£ 

ai 1 s 


®21 = 


■ ax ^ 

"3 

R* 

a p ^ 

G 


!ll 

ar 

» 

= 

apg 

r 

ap 

£2 . 

-< 

Tsr 

ap^ 

2G ‘ 

^ - 


2 r 

8 P 1 | 

r 

2 Pi 

ax 2 _ 

P 3 R 

• 

3Pl 

G 



^22 = T?: = 


^2 3r 

* - 

r apg 


g" ^3 '’2 


n» ‘’3^’ 

R - ^ sin p. 


^ . ^^2 ^2 ar 

= aTi; = r'aTi; 


'3 gTg + h: 



3 X 3 

^32 = TT- = 02 X 3 


H ^3 3r 

■« ■■■ 


G{G + H ) r a og 
19 





20 



and 


The derivatives 



are determined from equation (39). 


ar 

= 0 

k s 1 and 3 




aPk 





(7 

a r 

ap 2 

r 

" P 

-2»2 Jf- '■ ( 

i 

Q cos V 1 

Z;P|j ^ ^ 
2r^ 

1 ] 

a r 

r 

-2U p* 

rZ, 

Zl 

\1 Mn n 

3P|j 

P 

(2ojj)3''2 \/ir 

’ VS**!} 

2M0(2P4)^^7J 

ar 

ack 

= 0, 

1 

k s 3 and 11 



; ■ 

a r 

M •• 

pfc 





ao ^ 

F 

^2 0 




\J 3jr^ 

r 

-2»2 I"- ■’ ( 

f 

ZlPjjOgV 


aog 

" P 

Q cos 0 y 

\ 

2m2q / 


The derivatives 

and 

®”k 

a(E - ♦) 
ap^ - 

are determined from equation (4l) 


a(E - «) _ 

rr+ - e2) 


30i 


Ai n T - x^ *(® sin ♦) 

(1 + y1 - e*^ + e cos ♦) ■ ■ 

k 


3(E - ») 
a p\f 


~ e sin « 


-r/P 


\/l - e^ aH cos ♦) 
•> +■ 


a 0. 


a 0, 


(14, - e^) 

- e sin ^ 


(102) 


y. (T""'-" — S' x\ *(e sin ♦) 

(1 + y1 - + e cos ♦) 


3 \/l - a(e cos ♦) 


8 Pt 


3 0 ,. 


k = 1, 2, 3, ^1 


. j , . a(e cos «) a(e sin ♦) a(e oos ^) . 

The derivatives rr > rz 1 rr 1 and 


a Ob 


a 0 ,. 


a Pi. 


a(e sin derived from equations ( 3 J)> and ( 32 ). 

,73‘>k V'V 


a(e :?in 


= 0 , k = 1 and 3 


= 0 sin a, - 

ap? 2v^0 


a(e ain *) 


2 2pij 2 u0(8p4)^''^ 


a(e 3in ») _ ^ ^ 3 a,vv4 

3 01, '■ -- 


a(e sin -i) 


= e cos $ 


= Q cos 01 - Z2 -V 
aop ^ 2v^Q 


- 0 , k = 1 and 3 

8pj, 


a(e cos 
aP2 


ZiP4P2 

= 0 cos 01 - ^ — 


a(e cos 
ao2l 


- h \ 


2PI| 2mO(8o4) 


a (e c6s ii> ) 


= 0 , k = 3 and M 


a ( e cos ■» ) 


e sin lb 


a(e cos 41 ) 


/ Z 1 P 402 

= - 0 sin oi + 


ajl ~ e‘ 


And, finally, the derivatives — and 


a til - e‘ 


are derived 











7.0 NUMERICAL APPLICATIONS 

!(' V 

In reference 2, a number of numerical experiments were performed, 
comparing the PS(u) formulation with the KS formulation described in ref- 
erence 4; - In. reference 2, a complete' description of the force models 
and the initial conditions for different orbits were Riven. These same 
experiments_were performed except that the orbits were integrated with 
the upe- of F3(*) elements. 


7.1 Example 1 

The first example is a highly eccentric orbit (e = 0.95) about the 
Earth. The satellite is subject to the perturbing conservative potential 
of the Earth's oblateness (Jg) and -the perturbing forces due to the Moon, 

As in reference 2 , the PS(<i>) elements were integrated with an RK-45 fixed- 
step method. The PS{u), KS, and PsC'l’} solutions, at the end of 50 revolu- 
tions, were compared to an extremely accurate reference solution found 
in reference 4. 

The problem description for the first example is as follows. 


Coordinate svstem: Xi X fixed in Earth equatorial plane. 

' 2 

Xg perpendicular to Earth equatorial plane. 
Initial conditions: X| = 0.0 km, X 2 = -5888.9727 km 

X 3 = -3400.0 km 

ky = 10.691338 km/sec 

Xg = 0 km/sec, X 3 = 0 km/sec 

The time of comparison is at 288.12768941 days. The Earth oblateness 
and lunar perturbation models are described in section 9 . 0 , parts a and 
b of reference 2 . 




TABLE I.- EXAMPLE 1 


xi , km xp» km xo, km Energy Steps/rev 

check 


Ref. 

-24 219.0503 

227 -962.1064 

129 753*4424 





PSu 

8;2157 

.0072 

.3843 

0.21 

X 

10“® 

40 

KS f 

8.2221 

.0079 

.3883 

,20 

X 

10"S 

40 

PS‘> 

43.7435 

99.6752 

74.1071 

0.33 

X 

10-7 

40 

PSu 

.0422 

.1066 

.4421 

.22 

X 

10 *® 

80 

KS 

.0422 

.1070 

.4420 

.11 

X 

Iq -8 

80 

PS'I> 

21.2402 

3.4532 

4.2147 

.98 

X 

10-9 

80 

PSu 

.0499 

..1064 

.4424 

.18 

X 

10-9 

160 

KS - 

.0499 

II 

.1064 

..4424 

.14 

X 

10-9 

160 

PS 'll 

.1736 

.1694“ 

.4783 

.42 

X 

10-10 

160 


Note: The values shown are identical to the corresponding reference 

value except for the digits shown in the table. The energy check is the 
difference between the integrated value for the total energy and the total 
energy computed from position, velocity, and perturbing potential, divided 
by the integrated energy. 

As one can see, the PS(u) and KS solutions produce similar results, 
whereas the PS( *) solutions produces results that are not as accurate 
for the same step size. The .reason for this is clear. When a highly 
eccentric satellite approaches the apogee point, the lunar perturbations 
become very pronounced. The PS(u) and KS both have independent variables 
that are similar to the eccentric anomaly, whereas the PS( *) has the true 
anomaly as the independent variable. Therefore, near the apogee, the 
PS(u) and KS formulations are evaluating the strong lunar perturbations 
more often than the PS( ♦) formulation for the simple reason of geometry. 


7.2 Example 2 

In the second example, a satellite is subject to the J 2 perturba- 
tions and the drag perturbations of the form described in reference 3 , 
section 9.0 a and b. The satellite is in orbit about the Earth, with 
an eccentricity of e = 0.174. Again, the experiments were performed 
in the same manner as described inn the reference. One can see from Table 
II that all formulations Integrate this example easily. 

The problem description'Tor the second example is as follows. 
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Coordinate SYstem: same as that for example 1 

Initial conditions: x-| = 0.0 km, xg = -5888.9727 km, X 3 r,= -3*100.0 km 

xy = 8.3 km/sec, xg = 0.0 km/sec, Xg = 0.0 km/sec 

The time of comparison is at 4.30899150 days, after approximately 50 revo- 
lutions. The oblateness and drag perturbations were computed from the 
equations in parts a and c of section 9, reference 2. 

TABLE II.- EXAMPLE 2 





x-| , km 

X 2 , km 

x^ , km 
3 

Energy 

Steps, 

Ref. 

819.9225 

-5960.6168 

-3175.3500 

0.70 X 

10-12 

160 

PSu 

.9225 

.6163 

.3500 

.27 X 

Iq-8 

30 

KS 

.9225 

.6168 

-.3500 

X 

10-® 


PS 4i 

.9225 

.6168 

.3500 

i 

.7^4 X 

10-5 

30 

PSU 

.9225 

.6168 

• 3500 

. 22 

10-9 

50 

KS 

.9225 

V .6168 

.3500 

. 18 X 

10-9 

50 

PS ♦ 

.9225 

.6168 

.3500 

.70 X 

10-10 

50 


7.3 Example 3 

In this last example, a particle is placed at the libration 
in the Earth-Moon system. If the Earth and the Moon are assumed 
in circular orbits about their center of mass and if an infinite 
mass is placed 60 ‘ ahead of or behind the Moon at the same radiu 
circular speed as the Moon, then from the theory of the restrict 
body problem, the mass will remain at either of these points ind 
If the oarticle is placed there with a small initial position er 
the particle will librate about that point. The motion of such 
is computed by the PS ( 'll) , PS(u), and KS methods. As before, the 
is the same as described in reference 2, section 9.3* 


point 
to move . 
simal 

3 and ; 

ed three- 
efinitely. 
ror, then 
a particle 
experiment 


The problem description is as follows. 

[/ I: >' 

Coordinate system; ; Xi and X 2 ,/;in Earth-Moon plane; X.^^, perpendicular 
Initial conditions: x^ = 192 3P0 km X 2 = 332 900.1655 km, x^ = 0 km 

x-| = -0.8872840638 km/sec, X2 = 0.5125402247 km/sec 
Xj = 0; km/sec 





The final time is 10 
is given in part d, 

000 days. The lunar ephemeris 
section 9, of reference 2. 

used 

in this 

example 



TABLE III 

.- EXAMPLE 3 






x-| , km 

X 2 » km 

r , km 

Energy 

check 

V 

Steps/rev 

Ref. 

-178 094.8956 

-340 462.3726 

384 229.6435 0 

.50 

X 


80 

PSu 

5.4693 

1.9627 

.5463 

• 12 

X 

lO**"^ 

10 

KS 

6,4119 

.8837 

30.7993 

.72 

X 

10-5 

10 

PS* 

__^?.0474 

.0933 

.4664 - 

.12 

X 

10-<^ 

10 

PSu 

.9156 

.3688 

.6495 

.25 

X 

10-9 

' 20 

KS 

.9470 

.3989 

.6907 - 

// 

.24 

X 

,0-6 

20 

PS* 

.9036 

.3753 

1' .6497 

.25 

X 

10-9 

20 

PSu 

.8963 

.3726 

)] .6439 

J/ 

.36 

X 

10-11 

40 

KS 

.8972 

.3735 / 

.6451 

.75 

X 

10-5 

40 

PS* 

.8958 


.6439 

.38 

X 

10-9 

40 


From table III, one finds that the PS(<^) and PS(u) formulations give more 
accurate solutions than the KS method, for step ^dzes of 10, 20, and 40 
steps /rev. 
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8/0“ CONCLUSIONS 


The differential equations for the perturbed two-body elliptical 
motion are defined in elements that ore similar to the two-body classical 
canonical Poincar^ elements, which have time as the independent variable. 
A new Hamiltonian is defined that enables the introduction of two more 
canonical variables, the total energy and the time element. Also, these 
new elements have as their independent variable, the generalized true 
anomaly. 

The accuracy of the solution of this new set is equal to that of 
the PS(u) or KS formulations except for the highly 'eccentric case where 
the PS( di) is at a disadvantage because of geometric reasons< 
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APPENDIX 


COMPUTATIONAL PROCEDURE 


In the numerical integration of a trajectory by means of the PS(*) 
elements, one may follow the procedure described Jelow. Initially, the 
position and velocity vectors are converted bo the PS(0 elements (COTOPS) 
The numerical integrator will evaluate the PS(*) differential equations 
(RSDEEO) and estimate the value of the PS(|^) elements at t s tq + at, 
the independent variable. In PSDEBO, one must evaluate the value of the 
perturbing forces that are usually given as a function of the Cartesian 
state vector and time. Therefore, PSDEEO must convert the PS* elements 
into the coordinates (PSTOCO). The algorithms for COTOPS, PSDEEO, and 
PSTOCO are outlined in the following sections. The left column gives 
the quantity to be comouted, and the right column references the equation 
number in the text. 


COTOPS 

Given X, x, and time, transform to £. and £. 
Evaluate the potential, V. 

Then sequentially compute: 

• = L 

Gi, Gg, G3 

G 3 

H 

^ P ^ ^ 


* ^ °3 S .. .. 

* P 3 

R ~ ' 

r cos 
r sin ''•) 

. , ^ “ -^ECEDflMG PAGE BL.AWK NOT FILMED 

r 


From equation 
(75) 

(70)-(72) 

(73) 

(7iO 

(77) 

(79) 

(80) 

; (81) 

(87) 

( 88 ) 

(89) 

. (90) 

( 91 ) 

(92) 

(85) 
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i 

II 


n 

Given Oj a, transform to x, x, b = time, 
Seauentially compute: 

z, 

Z2 

0 

e cos ^1 

e sin 

P 

r 

R 

G 

H = G - (G - H) 

xi. X2, X3 

\j] - 

E - * 

t 

a 

p 

R» 

Xi, Xg, X3 


Zl 

Z 2 

E - ♦ 

G. 

PSTOCO 


( 82 ) 

(8H) 

(^3) 

(9»l) 

( 96 ) 

( 97 ) 
( 95 ) 


From equation 
( 33 ) 

( 3 *») 

( 35 ) 

(31) 

( 32 ) 

( 26 ) 

( 39 ) 

( 38 ) 

(83) 

(2i») 

■ (36) 

( 20 ) 

(M) 

(ilO) 

( 9 ) 

( 6 »») 

( 69 ) 

( 66 )- ( 68 ) 
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Given n, conjDUfce 


da 

dx' 


da 

dt* 


Proceed to PSTOCO, 



